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ABSTRACT 

Equilibrium models of differentially rotating nascent neutron stars are constructed, 
which represent the result of the accretion induced collapse of rapidly rotating white 
dwarfs. The models are built in a two-step procedure: (1) a rapidly rotating pre- 
collapse white dwarf model is constructed; (2) a stationary axisymmetric neutron star 
having the same total mass and angular momentum distribution as the white dwarf is 
constructed. The resulting collapsed objects consist of a high density central core of 
size roughly 20 km, surrounded by a massive accretion torus extending over 1000 km 
from the rotation axis. The ratio of the rotational kinetic energy to the gravitational 
potential energy of these neutron stars ranges from 0.13 to 0.26, suggesting that some 
of these objects may have a non-axisymmetric dynamical instability that could emit 
a significant amount of gravitational radiation. 
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1 INTRODUCTION 

The accretion induced collapse of a rapidly rotating white 
dwarf can result in the formation of a rapidly and dif- 
ferentially rotating compact object. It has been suggested 
that such rapidly rotating objects c ould emit a s ubstan- 
tial amount of gravitational radiation (Thorne 1995), which 



might be observable by the gravitational wave observatories 
such as LIGO, VIRGO and GEO. It has been demonstrated 
that if the collapse is axisymmetric, the energy emitted by 
gravi tational waves is rat her small ( Miiller & Hillebrandt 



198l|; [Finn & Evans 1990] [Monchmeyer ct al. 199 If |Zwcrgci 



& Miiller 1997). However, if the collapsed object rotates 



than the dynamical timescale of the system. An interest- 
ing class of secular and dynamical instabilities only occur 
in rapidly rotating stars. One convenient measure of the ro- 
tation of a star is the parameter /3 — T Iot /\W\, where T rot 
is the rotational kinetic energy and W is the gravitational 
potential energy. Dynamical and secular instabilities set in 
when f3 exceeds the critical values fid and (3 S respectively. 
It is well known that Pd ~ 0.27 and j3 a ~ 0.14 for uni- 
formly rotating, constant density and incompr essible stars, 
the Maclaurin spheroids ( |Chandrasekhar 196!: ). Numerous 
numerical simulations in Newtonian theory show that (3d and 
p a have roughly these same values for differentially rotating 
polytropes with the same specific a ngular momentum distri- 
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axisymmetric perturbations of the form e im * , where ^ is for Polytropes with some other angular momentum distri- 
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is drive n by hydrodynamics and gravity , and develops on a 
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processes such as viscosity or gravitational radiation reac- 
tion, and the growth time is determined by the dissipative 
timescale. These secular timescales are usually much longer 



Most of the stability analyses to date have been carried 
out on stars having simple ad hoc rotation laws. It is not 
clear whether these rotation laws are appropriate for the 
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nascent neutron stars formed from the accretion induced 
collapse of rotating white dwarfs. 

New-born neutron stars resulting from the core col- 
lapse of massive stars with realistic rotation laws were stud- 
ied by Monchmeyer, Janka and Miiller (Monchmeyer & 
Miiller 1988; Jank a & Monchmeyer 1989ab), and Zwerger 
and Miiller (1997). The study of Monchmeyer et al. shows 
that the resulting neutron stars have (3 < 0.14. Zwerger and 
Miiller carried out simulations of 78 models using simpli- 
fied analytical equations of state (EOS). They found only 
one model having /3 > 0.27 near core bounce. However, 
f3 remains larger than 0.27 for only about one millisecond, 
because the core re-expands after bounce and slows down. 
The pre-collapse core of that model is the most extreme one 
in their large sample: it is the most rapidly and most dif- 
ferentially rotating model, and it has the softest EOS. In 
addition, they found three model s with 0.14 < [3 < 0.27. 
Rampp, Miiller and Ruffert (1998) subsequently performed 
3D simulations of three of these models. They found that 
the model with /3 > 0.27 shows a non-linear growth of a 
non-axisymmetric dynamical instability dominated by the 
bar mode (m = 2). However, no instability is observed for 
the other two models during their simulated time interval 
of tens of milliseconds, suggesting that they are dynami- 
cally stable. Their analysis does not rule out the possibility 
that these models have non-axisymmetric secular instabili- 
ties, because the secular timescale is expected to range from 
hundreds of milliseconds to few minutes, much longer than 
their simulation time. 

The aim of this paper is to improve Zwerger and 
Miiller's study by using realistic EOS for both the pre- 
collapse white dwarfs and the collapsed stars. For the pre- 
collapse white dwarfs, we use the EOS of a zero-temperature 
degenerate electron gas with electrostatic corrections. A hot, 
lepton rich protoneutron star is formed as a result of the 
collapse. This protoneutron star cools down to a cold neu- 
tron star in about 20 s (see e.g. Burrows & Lattimer 1986), 
which is much longer than the dynamical timescle. So we 
adopt two EOS for the collapsed stars: one is suitable for 
protoneutron stars; the other is one of the standard cold 
neutron-star EOS. 

Instead of performing the complicated hydrodynamic 
simulations, however, we adopt a much simpler method. 
We assume (1) the collapsed stars are in rotational equi- 
librium with no meridional circulation, (2) any ejected ma- 
terial during the collapse carries a negligible amount of mass 
and angular momentum, and (3) the neutron stars have the 
same specific angular momentum distributions as those of 
the pre-collapse white dwarfs. The justifications of these as- 
sumptions will be discussed in Section ^[ Our strategy is as 
follows. First we build the equilibrium pre-collapse rotating 
white dwarf models and calculate their specific angular mo- 
mentum distributions. Then we construct the resulting col- 
lapsed stars having the same masses, total angular momenta 
and specific angular momentum distributions as those of the 
pre-collapse white dwarfs. All computations in this paper 
are purely Newtonian. In the real situation, if a dynamical 
instability occurs, the star will never achieve equilibrium. 
However, our study here can still give a useful clue to the 
instability issue. 

The paper is organized as follows. In the next section 
we present equilibrium models of pre-collapse, rapidly and 



rigidly rotating white dwarfs. In Section |§| we construct the 
equilibrium models corresponding to the collapse of these 
white dwarfs. The stabilities of the collapsed objects are dis- 
cussed in Section ^. Finally, we summarize our conclusions 
in Section |E| 



2 PRE-COLLAPSE WHITE DWARF MODELS 
2.1 Collapse mechanism 

As mass is accreted onto a white dwarf, the matter in the 
white dwarf's interior is compressed to higher densities. 
Compression releases g ravitational en ergy and some of the 
energy goes into heat ( Nomoto 1982[ ). If the accretion rate 
is high enough, the rate of heat generated by this compres- 
sional heating is greater than the cooling rate and the cen- 
tral temperature of the accreting white dwarf increases with 
time. 

The inner core of a carbon-oxygen (C-O) white dwarf 
becomes unstable when the central density or temperature 
becomes sufficiently high to ignite explosive carbon burn- 
ing. Carbon deflagration releases nuclear energy and causes 
the pressure to increase. However, electron capture behind 
the carbon deflagration front reduces the temperature and 
pressure and triggers collapse. Such a white dwarf will ei- 
ther explode as a Type la supernova or collapse to a neu- 
tron star. Which path the white dwarf takes depends on the 
competit ion between th e nuclear energy release and electron 
capture (Nomoto 1987). If the density at which carbon ig- 



nites is higher than a critic al density of about 9 x 10 9 g cm 3 
([Timmes fc Woosley 1992), electron capture takes over and 



the white dwarf will collapse to a neutron star. However, if 
the ignition density is lower than the critical density, carbon 
deflagration will lead to a total disruption of the whole star, 
leaving no remn ant a t all. More recent calculations by Bravo 
& Garcia-Senz ( |l999| ), taking into account the Coulomb cor- 
rections to the EOS, suggests that this critical density is 
somewhat lower: 6 x 10 9 g cm -3 . The density at which car- 
bon ignites depends on the central temperature. The cen- 
tral temperature is determined by the balance between the 
compressional heating and the cooling and so strongly de- 
pends on the accretion rate and accretion time. For zero- 
temperature C-0 white dwarfs, carbon ignites at a density 
of about 10 10 g cm' 3 (ISalpctcr fe Van Horn 1969| |Ogata' 



Iyetomi & Ichimaru 1991), which is higher than the above 



critical density. If the accreting white dwarf can somehow 
maintain a low central temperature during the whole accre- 
tion process, carbon will ignite at a density higher than the 
critical density, and the white dwarf will collapse to a neu- 
tron star. The fate of an accreting white dwarf as a function 
of the accretion rate and the white dwarf's initial mass is 



summarized in two diagrams in the paper of Nomoto (1987 



(see also Nomoto & Kondo 1991). Roughly speaking, low ac- 
cretion rates (M ^ 1O _8 M0 yr _1 ) and high initial mass of 
the white dwarf (M <; 1.1 Mq), or very high accretion rates 
(near the Eddington limit) lead to collapse rather than ex- 
plosion. 

Under certain conditions, an accreting oxygen-neon- 
magnesium ( O-Ne-Mg) white dwarf can also collaps e to a 



neutron star (Nomoto 1987; Nomoto fc Kondo 1991). The 



collapse is triggered by the electron captures of ^Mg and 
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20 Ne at a density of 4 x 10 9 g cm" 3 . Electron captures not 
only soften the EOS and induce collapse, but also gener- 
ate heat by 7-ray emission. When the star is coll apsed to 



a centr al densi ty of 10 g cm , oxygen ignites (Nomoto 



& Kondo 1991). At such a high density, however, electron 
captures occur at a faster rate than the oxygen burning and 
the white dwarf collapses all the way to a neutron star. 

In this section, we explore a range of possible pre- 
collapse white dwarf models. We assume that the white 
dwarfs are rigidly rotating. This is justified by the fact that 
the timescale for a magnetic field to suppress any differential 
rotation, tb, is short compared with the accretion timescale. 
For example, tb ~ 10 3 years if the massive white dwarf has 
a magnetic field B ~ 100G. We construct three white dwarf 
models using the EOS of a zero-temperature degenerate 
electro n gas with Coulomb corrections derived by Salpeter 
( 1961 ). All three white dwarfs are rigidly rotating at the 
maximum possible angular velocities. Model I represents a 
C-0 white dwarf with a central density of p c = 10 10 g cm -3 , 
the highest p c a C-0 white dwarf can have before carbon ig- 
nition induces collapse. Model II is also a C-O white dwarf 
but has a lower central density, p c = 6 x 10 9 g cm -3 . This is 
the lowest central density for which a white dwarf can still 
collapse to a neutron star after carbon ignition. Model III 
is an O-Ne-Mg white dwarf with p c — 4 x 10 9 g cm" 3 , that 
is the density at which electron captures occur and induce 
the collapse. Since the densities are very high, the pressure 
is dominated by the ideal degenerate Fermi gas with elec- 
tron fraction Z/A = 1/2 that is suitable for both C-O and 
O-Ne-Mg white dwarfs. Coulomb corrections, which depend 
on the white dwarf composition through the atomic number 
Z, contribute only a few per cent to the EOS at high den- 
sities, so the three white dwarfs are basically described by 
the same EOS. 



2.2 Numerical method 

We treat the equilibrium rotating white dwarfs as rigidly 
rotating, axisymmetric, and having no internal motion other 
than the motion due to rotation. The Lichtenstein theorem 
(see Tassoul 1978) guarantees that a rigidly rotating star 
has reflection symmetry about the equatorial plane. We also 
neglect viscosity and assume Newtonian gravity. Under these 
assumptions the equilibrium configuration is described by 
the static Euler equation: 



v ■ Vw 



VP 

P 



(1) 



where P is pressure; p is density; $ is the gravitational po- 
tential, which satisfies the Poisson equation 



V 2 $ = AivGp , 



(2) 



where G is the gravitational constant. The fluid's velocity 
v is related to the rotational angular frequency f2 by v — 
Qzue^, where vj is the distance from the rotation axis and 
ecp is the unit vector along the azimuthal direction. The EOS 
we use is barotropic, i.e. P = P(p), so the Euler equation 
can be integrated to give 



where C is a constant. The enthalpy (per mass) h is given 
by 



h= — 



dP 
P 



(4) 



and is defined only inside the star. The boundary of the star 
is the surface with h = 0. 

The equilibrium configuration is determin ed by 



Hachisu's self-consistent field method (Hachisu 1986): given 
an enthalpy distribution hi, we calculate the density distri- 
bution pi from the inverse of equation (^) and from the EOS. 
Next we calculate the gravitationalpotential $i everywhere 
by solving the Poisson equation (fe|). Then the enthalpy is 
updated by 



hi+i — Cia 



(•») 



with Ci+i and flf+j det ermin ed by two boundary condi- 



tions. In Hachisu's paper (1986), the axis ratio, i.e. the ratio 
of polar to equatorial radii, and the maximum density are 
fixed to determine Cj+i and fi 2 +1 . However, we find it more 



convenient in our case to fix the equatorial radius R e 
central enthalpy h c , so that 



C, 



Or 



-—[a +1 - <f t (A)} , 

tie 



and 

(6) 
(7) 



where &i(0) and &i(A) are the gravitational potential at 
the centre and at the star's equatorial surface respectively. 
The procedure is repeated until the enthalpy and density 
distribution converge to the desired degree of accuracy. 

We used a spherical grid with L radial spokes and 
N evenly spaced grid points along each radial spoke. The 
spokes are located at angles 9i in such a way that cos 0i cor- 
respond to the zeros of the Legendre polynomial of order 
2L — 1: P2L-1 (cos #i) = 0. Because of the reflection sym- 
metry, we only need to consider spokes lying in the first 
quadrant. Poisson's equation is s olved using the technique 
described by Ipser and Lindblom (199C). The special choice 
of the angular positions of the radial spokes and the fi- 
nite difference scheme make our numerical solution equiva- 
lent to an ex pansion in Legendre poly nomials through order 
/ = 2L - 2 ( |lpscr fc Lindblom 1990| ). Although the white 
dwarfs we consider here are rapidly rotating, the equilib- 
rium configurations are close to spherical, as demonstrated 
in the next subsection. So a relatively small number of ra- 
dial spokes are adequate to describe the stellar models accu- 
rately. We compared the results of (L, N) = (10, 3000) with 
(L, N) — (20, 5000) and find agreement to an accuracy of 
10" ° . The accuracy of the model can also be measured by 
the virial theorem, which states that 2T ro t + W + 311 = for 
any equilibrium star (see e.g. Tassoul 1978). Here T rot is the 
rotational kinetic energy; W is the gravitational potential 



energy and Tl = J Pd 3 x. We define 



2T rot + W + 3H 



W 



All models constructed in this section have e ~ 10 



(8) 



h = C ■ 



(3) 
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Figure 1. Meridional density contours of the rotating white 
dwarf of Model I. The contours, from inward to outward, cor- 
respond to densities p/p c =0.8, 0.6, 0.4, 0.2, 0.1, 10~ 2 , 10 -3 , 
10" 4 , 10~ 5 and zero. 

2.3 Results 

We constructed three models of rigidly rotating white 
dwarfs. All of them are maximally rotating: material at 
the star's equatorial surface rotates at the local orbital fre- 
quency. Models I and II are C-O white dwarfs with cen- 
tral densities p c — 10 10 g cm -3 and p c — 6 x 10 9 g cm~ 3 
respectively; Model III is an O-Ne-Mg white dwarf with 
p c = 4x 10 9 g cm -3 . The properties of these white dwarfs are 
summarized in Table |l|. We see that the angular momentum 
J decreases as the central density p c increases, because the 
white dwarf becomes smaller and more centrally condensed. 
Although T rot increases with p c , \W\ increases at a faster 
rate so that f3 = T rot /|jy| decreases with increasing p c . We 
also notice that the mass does not change much with increas- 
ing p c . The reason is that massive white dwarfs are centrally 
condensed so their masses are determined primarily by the 
high density central core. Here the degenerate electron gas 
becomes highly relativistic and the Coulomb effects are neg- 
ligible, so the composition difference is irrelevant. Hence the 
white dwarf behaves like an n — 3 polytrope, whose mass in 
the non-rotating case is independent of the central density. 

The masses of our three models are all greater than 
the Chandrasekhar limit for non-rotating white dwarfs. A 
non-rotating C-O white dwarf with p c = 10 10 g cm -3 has a 
radius R = 1300 km and a mass M — 1.4OM0. When this 
white dwarf is spun up to maximum rotation while keeping 
its mass fixed, the star puffs up to an oblate figure of equa- 
torial radius R e — 4100 km and polar radius R p — 2700 km, 
and its central density drops to p c = 5.5 x 10 s g cm -3 . This 
peculiar behaviour is caused by the soft EOS of relativistic 
degenerate electrons, which makes the star highly compress- 
ible and also highly expansible. When the angular velocity 
of the star is increased, the centrifugal force causes a large 
reduction in central density, resulting in a dramatic increase 
in the overall size of the star. 

Figures [l]-^ display the density contours of our three 
models. The contours in the high density region remain more 
or less spherical even though our models represent the most 
rapidly rotating cases. The effect of rotation is only to flatten 
the density contours of the outer region in which the density 
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Figure 2. Same as Figure Fl but for Model II 



z/km 




1000 2000 

m/km 



Figure 3. Same as Figure |l| but for Model III 

is relatively low. This suggests that the white dwarfs are 
centrally condensed, and is clearly demonstrated in Figure ^, 
where the cylindrical mass fraction 

rrivj = — / dvo' vj I dz' p(vj , z) (9) 

JO J-oo 

is plotted. In all of our three models, more than half of the 
mass is concentrated inside the cylinder with vj « 0.2R e . 

Figure shows the specific angular momentum j as a 
function of the cylindrical mass fraction m m , normalized so 
that J j(m-ro) drtito = 1. The j(m II ,)-curves for the three 
models are almost indistinguishable except in the region 
where m ra ~ 1. The spike of the curve near m m — 1 can be 
understood from Figure ^, where we see that w 1 when 
vjjR e <; 0.6. However, j = (M/J)Qzu 2 oc -cu 2 . These two 
make the values of j in the interval 0.6 2 ^ j/jirn^ = 1) < 1 
squeeze to the region m w 1, and the spike results. We shall 
point out in the next section that this spike causes a seri- 
ous numerical problem in the construction of the equilib- 
rium models of the collapsed objects. The problem can be 
solved by truncating the upper part of the j(m m ) curve. 
The physical justification is that the material in the outer 
region contributes only a very small fraction of the total 
mass and angular momentum of the star, as illustrated in 
Table El for Model I. The situations for the other two models 
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Table 1. The central density p c , mass M, angular momentum J, rotational frequency Q, rotational kinetic energy T ro t, the ratio of 
rotational kinetic to gravitational energies /3, equatorial radius R e and polar radius R p of three rigidly and maximally rotating white 
dwarfs. 
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Figure 4. Cylindrical mass fraction function of to. Solid 

line corresponds to Model I; dotted line, Model II, and dashed 
line, Model III. 
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Figure 5. Normalized specific angular momentum j as a function 
of the cylindrical mass fraction m ro . The curves for the three 
models are indistinguishable except in the region very close to 
m m = 1, which is magnified in the inset. 



are very similar and so are not shown. We see that material 
in the region where vj j R e > 0.9 [j(m m ) / > 0.81] con- 
tributes less than 10~ 5 of the total mass and 10" 4 of the 
total angular momentum. So the upper 19 per cent of the 
j(m-zj)-curve has little influence to the inner structure of the 
collapsed star. While this region is important for the struc- 
ture of the star's outer layers, that part of the star is not of 
our primary interest since the mass there is too small to de- 



Table 2. The outer layers of Model I white dwarf. is the 
angular momentum of the material inside the cylinder of radius 
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velop any instability that can produce a substantial amount 
of gravitational radiation. 



3 COLLAPSED OBJECTS 

In this section, we present the equilibrium new-born 
neutron-star models that may result from the collapse of 
the three white dwarfs computed in the previous section. 
Instead of performing hydrodynamic simulations, we adopt 
a simpler approach: 

First, we assume the collapsed stars are axisymmetric 
and are in rotational equilibrium with no meridional circu- 
lation. Second, we assume the EOS is barotropic, P = P(p). 
These two assumptions imply that (1) the angular veloc- 
ity Q is a function of zu only, i.e. dSl/dz = 0, and (2) the 
Solberg condition is satisfied, which states that dj/dzu > 
for stable barotropic stars in rotational equilibrium (see e.g. 
Tassoul 1978). The angular velocity profile (dQ/dz = 0) 
is observed in the simulations of Monchmeyer, Janka and 
Miiller (Monchmeyer & Miiller 1988; Janka & Monchmeyer 
1989ab). Third, we are only interested in the structure of 
the neutron stars within a few minutes after core bounce. 
The timescale is much shorter than any of the viscuous 
timescales, so viscosity does not have time to change the 



angu l ar momentum of a fluid particle (Sutler fc Lindblom 



1987; sawyer 1989 



Lindblom & Detweiler 19791; van den 



Horn fc van Weert 1981| poodwin fc Pethick 198^ ). Finally, 



we assume no material is ejected during the collapse. It fol- 
lows, from the conservation of j and the fact that j is a 
function of vj only before and after collapse, that all par- 
ticles initially located on a cylindrical surface of radius vji 
from the rotation axis will end up being on a new cylindri- 
cal surface of radius zo2- And the Solberg condition ensures 
that all particles initially inside the cylinder of radius vj\ will 
collapse to the region inside the new cylinder of radius W2- 
Hence the specific angular momentum distribution j(m^) of 
the new equilibrium configuration is the same as that of the 
pre-collapse white dwarf; here m ro is the cylindrical mass 
fraction defined by equation (H). 
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Based on these assumptions, we constructed equilib- 
rium models of the collapsed objects with the same masses, 
total angular momenta and ]'(mw) as the pre-collapse white 
dwarfs. 



3.1 Equations of state 

The gravitational collapse of a massive white dwarf is halted 
when the central density reaches nuclear density where the 
EOS becomes stiff. The core bounces back and within a few 
milliseconds, a hot (T <; 20 Mev), lepton rich protoneu- 
tron star settles into hydrodynamic equilibrium. During the 
so-called Kelvin-Helmholtz cooling phase, the temperature 
and lepton number decrease due to neutrino emission and 
the protoneutron star cools to a cold neutron star with 
temperature T < 1 Mev after several minutes. Since the 
cooling timescale is much longer than the hydrodynamical 
timescale, the protoneutron star can be regarded as in quasi- 
equilibrium. 

The EOS of a protoneutron star is expressed in the form 
P = P(p; s, Ye), where s and Y e are the entropy per baryon 
and lepton fractio n res pectively. As pointed out by Strobel, 
Scraab & Weigel ( 1999), the structure of a protoneutron star 
can be approximated by a constant s and Y e throughout the 
star, resulting in an effectively barotropic EOS. 

We used two different EOS for densities above 
The first is one of the standard EOS for cold 



10 10 g cm -3 

neutron stars. We adopt the Bethe- Johnson EOS (Bethe & 
Johnson 1974) for densities above 10 14 g cm" 3 , and BBP 



EOS (iBaym, Bethe & Pethick 1971) for densities in the re- 

TH 



gion 10 11 g cm 3 



10 i4 g cm- 



ties of these collapsed stars are lower than 4 x 10 14 



It turns out that the densi- 
g cm" 

and ideas about the EOS in this range have not changed 
very much since 1970's . The second is the EOS LPNSy L04 
of Strobel et al (1999)Q This corresponds to a protoneu- 
tron star 0.5 — Is after core bounce. It has an entropy per 
baryon s = 2/cb and a lepton fraction Y e = 0.4, where ks 
is Boltzmann's constant. We join both EOS to that of the 
pre-collapse white dwarf for densities below 10 10 g cm -3 . 
Hereafter, we shall call the first EOS the cold EOS, and the 
second one, the hot EOS. 



3.2 Numerical method 

We compute the equilibrium structure by Hachisu's self- 
cons istent field method modifi ed so that j(m ro ) can be speci- 
fied (Smith & Centrella 1992). The iteration scheme is based 
on the integrated static Euler equation (|l|) written in the 
form 



(10) 



where C is the integration constant, and M and J are the 
total mass and angular momentum of the star respectively. 
Given an enthalpy distribution hi everywhere, the density 
distribution pi is calculated by the EOS and the inverse of 



* The tabulated EOS is obtained from 


httn: / / www.physik.uni- 


muenchcn.de / sektion / suessmann / astro / eos / . 



equation (|4j). Next we compute the mass Mi and cylindrical 
mass fraction m TOi i by 



Mi = 4ir I dzu'zu' I dz pi(vj' , z') 
Jo _ Jo 

Wlra.i = TT I dvj' Vj' J dz'pi{zu',z') 



M, 







(11) 



I) 



and solve the Poisson equation V 2 $i = AivGpi to obtain the 
gravitational potential We then update the enthalpy by 
equation (|ic|): 



h i+1 (zu, z) 



= C l+ i - $i{zu,z) 
J x 2 



M, 



1+1 



dzu 



(12) 



o 



with the parameters d+i and (Ji+i/Mi+i) 2 determined by 
specifying the central density p c and equatorial radius R e . 
The procedure is repeated until the enthalpy and density 
distribution converge to the desired degree of accuracy. 

To construct the equilibrium configuration with the 
same total mass and angular momentum as a pre-collapse 
white dwarf, we first compute a model of a non-rotating 
spherical neutron star, use its enthalpy distribution as an 
initial guess for the iteration scheme described above and 
build a configuration with slightly different p c or R e . Then 
the parameters p c and R e are adjusted until we end up with 
a configuration having the correct total mass and angular 
momentum. 

Two numerical problems were encountered in this pro- 
cedure. The first problem is that when the angular momen- 
tum J is increased, the star becomes flattened, and the it- 
eration often oscillates among two or more states without 
converging. This problem can be solved by using a revised 
i terati on scheme suggested by Pickett, Durisen & Davis 
( 1996| ), in which only a fraction of the revised enthalpy h i+ i, 
h' i+1 = (1 — 8)hi+i + Shi, is used for the next iteration. Here 
S < 1 is a parameter controlling the change of enthalpy. We 
need to use S > 0.95 for very flattened configurations, and 
it takes 100 — 200 iterations for the enthalpy and density 
distributions to converge. 

The second problem has to do with the spike of the 
j(m ro )-curve near m ro = 1 (see Figure ||). The slope is so 
steep t hat it makes the iteration unstable. As discussed in 
Section 



2.3 



the material in the region very close to m ro = 1 
contains a very small amount of mass and angular momen- 
tum, so we can truncate the last part of the j(m ro )-curve 
without introducing much error. Specifically, we set a pa- 
rameter j c < j(m^ — 1), compute a quantity m c which 
satisfies j(m c ) = j c . Then we use the specific angular mo- 
mentum distribution j(m m ) = j(rn-^ m c ) instead of j(m-&). 
Typically, we choose j c /j{l) = 0.81 so that 1 — m c ~ 10" 5 
(see Table ^). Hence the distributions j(m ro ) and j(m^) 
are basically the same except in the star's outermost region, 
which is unimportant to the inner structure the star, and 
presumably also unimportant for the star's dynamical and 
secular stabilities. We also tried several different values of 
j c and found that the change of physical properties of the 
collapsed objects (e.g. the quantities in Table j3|) are within 
the error due to our finite-size grid. Thus the truncation is 
also justified numerically. 

We evaluate these stellar models on a cylindrical grid. 
This allows us to compute the integrals in equations ( |ll"l ) 
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Table 3. The central density p c , radius of gyration R g , charac- 
teristic radius R* and ratio of rotational kinetic energy to gravi- 
tational energy (3 of the collapsed objects with the cold and the 
hot EOS. 





Pc 


R g 


R* 







g cm -3 


km 


km 




Model I (cold EOS) 


3.7 x 10 14 


63 


670 


0.230 


Model I (hot EOS) 


1.4 x 10 14 


67 


650 


0.139 


Model II (cold EOS) 


3.5 x 10 14 


78 


800 


0.246 


Model II (hot EOS) 


0.79 x 10 14 


85 


800 


0.137 


Model III (cold EOS) 


3.2 x 10 14 


94 


940 


0.261 


Model III (hot EOS) 


0.27 x 10 14 


110 


940 


0.127 



and ( p^ ) easily. We find it more convenient however, to solve 
the the Poisson equation for the gravitational potential on a 
spheri cal gr id using the method described by Ipser and Lind- 
blom (199C). We have verified that the potential obtained in 



this way agrees with the result obtained with a cylindrical 
multi-grid solver to within 0.5%. However, the spherical grid 
solver (including the needed transformation from one grid to 
the other) is much faster than the cylindrical grid solver. The 
accuracy of our final equilibrium models can also be mea- 
sured by the quantity e defined in equation (^). The values 
of e for models computed in this section are few times 10 -4 . 



3.3 Results 

Table ^| shows some properties of the collapsed objects re- 
sulting from the collapse of the three white dwarfs in Sec- 
tion H. We define the radius of gyration, R g , and the char- 
acteristic radius, i?*, of the star by 



MRt 



2 ,3 

pzu a x 



m^{zu^R,) = 0.999 



(13) 
(14) 



We see that R g and R* that result from the same initial 
white dwarfs are insensitive to the neutron-star EOS, while 
there is a dramatic difference in the central density p c and 
the ratio of rotational kinetic energy to gravitational poten- 
tial energy (3. The collapsed stars with the hot EOS have 
smaller p c and /3 than those with the cold EOS. In fact, 
the central densities of these hot stars are less than nu- 
clear density. It is well-known that a non-rotating star can- 
not have a central density in the sub-nuclear density regime 
(4 x 10 11 g cm" 3 <; p <: 2 x 10 14 g cm -3 ) because the EOS is 
too soft to render the star stable against gravitational col- 
lapse. It has been suggested that if rotation is taken into 
account, a star with a central density in this regime can ex- 
ist. Such stars are termed 'fizzlcrs' in the literature ( Shapiro 
k Lightman 19761 |Tohline 198^ |Eriguchi fc Muller 19851 



Miiller fc Eriguchi 1985 j ; Hayashi, Eriguchi fc Hashimoto 



199S ; jlmamura, Duriscn & Pickett 2000| ). However, these so- 
called fizzlers in our case can exist for only about 20 s before 
evolving to rotating cold neutron stars. In order to build a 
stable cold model in the sub-nuclear density regime, the col- 
lapsed star has to rotate much faster, which is impossible 
unless the pre-collapse white dwarf is highly differentially 
rotating. 

We mention in Section |l| that Zwerger and Muller (1997 ) 
performed 2D hydrodynamic simulations of axisymmetric 



rotational core collapse. Their pre-collapse models are ro- 
tating stars with n = 3 polytropic EOS, which is close to 
the real EOS of a massive white dwarf. All of their pre- 
collapse models have a central density of 10 10 g cm -3 (see 
their Table 1). The model A1B3 in their paper is the fastest 
(almost) rigidly rotating star, but its total angular momen- 
tum J and /3 are respectively 22 per cent and 40 per cent less 
than those of our Model I of the pre-collapse white dwarf, 
though both have the same central density. This suggests 
that the structure of a massive white dwarf is senstive to 
the EOS. Zwerger and Muller state in their paper that no 
equilibrium configuration exists that has /3 > 0.01 for the 
(almost) rigidly rotating case. This assertion is confirmed 
by our numerical code. Zwerger and Muller adopt a sim- 
plified analytical EOS for the collapsing core. At the end 
of their simulations, the models A1B3G1-A1B3G5, corre- 
sponding to the collapsed models of A1B3, have values of 
(3 less than 0.07, far smaller than the f3's of our collapsed 
model I (see Table ^|), indicating that the EOS of the col- 
lapsed objects also play an important role on the final equi- 
librium configurations (or that their analysis violates one of 
our assumptions). 

Figures ^|-^ show the density contours of the collapsed 
objects. We see that the contours of the dense central region 
look like the contours of a typical rotating star. As we go out 
to the low density region, the shapes of the contours become 
more and more disk-like. Eventually, the contours turn into 
torus- like shapes for densities lower than I0~ 4 p c - In all cases, 
the objects contain two regions: a dense central core of size 
about 20 km and a low density torus-like envelope extending 
out to 1000 km from the rotation axis. Since we truncate 
the j(m ro )-curve, we cannot determine accurately the actual 
boundary of the stars. The contours shown in these figures 
have been checked to move less than one per cent as the 
cutoff j c /j(l) is changed from 0.7 to 0.9. This little change 
is hardly visible at the displayed scales. 

Figure ^| shows the rotational frequency / = Q/2n as a 
function of zu, the distance from the rotation axis. We see 
that the cores of the cold models are close to rigid rota- 
tion. The rotation periods of the cores of the cold neutron 
stars are all about 1.4 ms, slightly less than the period of 
the fastest observed millisecond pulsar (1.56 ms). A further 
analysis reveals that / oc zu~ a in the region zu <; 100 km, 
where a ~ 1.5 for the cold models and a ~ 1.4 for the hot 
models. 

To gain an insight into the structure of the envelope, 
we define the Kepler frequency fix at a given point on the 
equator as the angular frequency required for a particle to be 
completely supported by centrifugal force, i.e. fix satisfies 
the equation Q^zu = g, where g is the ma gnit ude of grav- 
itational acceleration at that point. Figure lid plots Q/Q,k 
as a function of zu along the equator. For the cold models, 
the curves increase from 0.5 at the centre to a maximum of 
about 0.95 at zu « 35 km, then decrease to a local mini- 
mum of about 0.8, and then gradually increase in the outer 
region. The curves of the hot models, on the other hand, 
increase monotonically from about 0.4 at the centre to over 
0.7 in the outer region. In all cases, centrifugal force plays 
an important role in the structure of the stars, especially in 
the low density region. 

Figure fLlI plots the cylindrical mass fraction 
function of zu. In all cases the cores contain most of the stars' 
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Model I (cold EOS) 



Model I (hot EOS) 



200 



z/km 



-200 




300 600 

toVkm 



200 



z/km 



-200 




300 

toVkm 




100 




100 



Figure 6. Meridional density contours of the neutron stars resulting from the collapse of Model I white dwarf. The left graphs correspond 
to the cold EOS, and the right graphs, the hot EOS. The contours in the upper graphs denote, from inward to outward, p/p c = 
10 — 3 , 10 — 4 , 10 — 5 , 10 — 6 , and 10 — 7 . The contours in the lower graphs denote, from inward to outward, p/p c =0.8, 0.6, 0.4, 0.2, 0.1, 
10" 2 , 10" 3 and 10~ 4 . 
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Figure 7. Same as Figure H but for Model II. 
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Model III (cold EOS) 
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Model III (hot EOS) 
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Figure 8. Same as Figure H but for Model III. 



mass. Material in the region tu ij 200 km occupies only a 
few per cent of the total mass, but it is massive enough that 
its self-gravity cannot be neglected in order to compute the 
structure of the envelope accurately. The envelope can be 
regarded as a massive, self-gravitating accretion torus. The 
same structure is also observed in the core collapse simula- 
tions of Janka & Monchmeyer (1989b) and Fryer & Heger 
(private communication with Fryer 



but different total masses and angular momenta. If the final 
model we end up with is unstable against axisymmetric per- 
turbations, there must be at least one model in the sequence 
such that 



dM 
dp c 







(17) 



which signals t he onset of instability ( Bisnovatsyi-Kogan & 



of zu, where T IC ,t{zu) and W{zu) are the rotational kinetic 
energy and gravitational potential energy inside the cylinder 
of radius zu, i.e. 



Fi gure T|~shows /3 ro = T Iot (zu)/\W(zu)\ as a function Blinnikov 1974| ). Here M is the total mass and p c is the cen 



7rot(w) = 2ir dva m (flw 



dz' p{zu' , 



(15) 



W(zu) — 2ty dvo'vo' I dz' p(zu' , z')<t(zu' , z) . (16) 
Jo Jo 

The values of (3 m approach /3 when w <; 40 km for the cold 
EOS models and when zu !> 100 km for the hot EOS mod- 
els. This suggests that material in the region zu J> 100 km 
contains negligible amount of kinetic energy, and any in- 
stability developed in this region could not produce strong 
gravitational waves. 



4 STABILITY OF THE COLLAPSED OBJECTS 

We first consider axisymmetric instabilities, i.e. axisymmet- 
ric collapse. This stability is verified when we construct 
the models. Recall that we start from the model of a non- 
rot at ing_sj}hj3ri£al_s^ 



an initial guess to build a sequence of rotating stellar mod- 
els with the same specific angular momentum distribution 



tral density. The partial derivative is evaluated by keeping 
the total angular momentum J and specific angular mo- 
mentum distribution j(m^) fixed. We have verified that all 
of our equilibrium models in the sequencies satisfy q > 0. 
Hence they are all stable against axisymmetric perturba- 
tions. 

We next consider non-axisymmetric instabilities. We 
have /3 = 0.23 - 0.26 for the cold EOS models and fi = 
0.13-0.14 for the hot EOS models (Table|). The hot models 
are probably dynamically stable but may be secularly unsta- 
ble. However, since they are evolving to cold neutron stars 
in about 20 s and their structures are continually changing 
on times comparable to the secular timescale, we shall not 
discuss secular instabilities of these hot models here. 

The values of /3 for the three cold neutron stars are 
slightly less than the traditional critical value for dynamical 
instability, /3d ~ 0.27. This critical value is based on simula- 
tions of differentially rotating polytropes having the j(m-^)- 
distribution of Maclaurin spheroids. However, recent sim- 
ulations demonstrate that differentially rotating polytropes 
having other j(m CT )-distributions can be dynamically unsta- 



ble for values of as low a s 0.14 (Pickett, Durisen & Davis 
199f; Centrella et al. 2000). The equilibrium configurations 



of some of those unstable stars also contain a low density 
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Figure 9. Rotational frequency / as a function to for the cold Figure 10. The ratio Q/Qk along the equator as a function of 

models (upper graph) and the hot models (lower graph). The to for the cold models (upper graph) and the hot models (lower 

inset in each graph shows / in linear scale in the central region. graph). 



accretion disk like structure in the stars' outer layers. This 
feature is very similar to the equilibrium structure of our 
models. Hence a more detailed study has to be carried out 
to determine whether the cold models are dynamically sta- 
ble. 

The subsequent evolution of a ba r-unstable object has 



been studied for the past 15 years f Dm-isen j ^juigold_& 
Toh line[T98^; |Williams fc Tohline 1988| |Houser fc Centrella 



199C ; Pickett, Durisen & Davis 1996; Smith, Houser & Cen 



trella 1 396^ | IIouacr 1998 |; Hew, Ccnti-clla fc Tohline 200 6} 



( Friedman fc Schutz 1975, 197S ) . In our case, R ~ 35 km (see 
Figure |l|), ft w 4000 rad s" 1 and (3 w 0.24, so r bar ~ 0.1 s. 
Gravitational waves may also driv e the r-mode instability 
(Lindblom, Owen fc Morsink 199£). The timescale is esti- 



mated by 



7.3 s 



10 14 g cm" 



<)"( 



4000 rad s" 



(19) 



Imamufa, Duricon fe Pickott 200q ; prown 20001). Itlp-found Mendell fc Owen 199SQ 

' ' ' ■ efirtina' n for t na innpr 9 



for the I — 2 r-mod e at low temperatures ( Lindblom 



that a bar-like structure develops in a dynamical timescale. 
However, it is still not certain whether the bar structure 
would be persistent, giving rise to a long-lived gravitational 
wave signal, or material would be shed from the ends of the 
bar after tens of rotation periods, leaving an axisymmetric, 
dynamically bar-stable central star. 

Even if the cold neutron stars are dynamically sta- 
ble, they are subject to various secular instabilities. The 
timescale of the gravitational-wave-driven bar-mode insta- 
bility can be estimated by 



where p is the average density. In- 
serting p for the inner 20-km cores of the cold stars, we have 
T r ~ 10 s 3> Tbar- The evolution of the bar-mode secular 
instability has only been studied in detail for the Maclaurin 
spheroids. These objects evolve through a sequence of de- 
formed non- axisymmetric configurations eventually to set- 
t le down as a more slowly ro t ating stable axisym metric star 



(Lindblom fe Detweiler 1977; Lai fc Shapiro 1995). It is gen 



erally expected that stars having more realistic EOS will 
behave similarly. 



Tim 



o.i.f-?-rf 

\35kmJ V 



4000 rad s" 



0.1 



(18) 
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Figure 11. Cylindrical mass fraction function of for 

the cold models (upper graph) and the hot models (lower graph). 



Figure 12. The quantity /3 ro as a function of ro for the cold 
models (upper graph) and the hot models (lower graph). 



5 CONCLUSIONS 

We have constructed equilibrium models of differentially ro- 
tating neutron stars which model the end products of the ac- 
cretion induced collapse of rapidly rotating white dwarfs. We 
considered three models for the pre-collapse white dwarfs. 
All of them are rigidly rotating at the maximum possible an- 
gular velocities. The white dwarfs are described by the EOS 
of degenerate electrons at zero t empe rature with Coulomb 
corrections derived by Salpeter (1961). 



We assumed that (1) the collapsed objects are axisym- 
metric and are in rotational equilibrium with no meridional 
circulation, (2) the EOS is barotropic, (3) viscosity can be 
neglected, and (4) any ejected material carries negligible 
amounts of mass and angular momentum. We then built the 
equilibrium models of the collapsed stars based on the fact 
that their final configurations must have the same masses, 
total angular momenta and specific angular momentum dis- 
tributions, j(m^), as the pre-collapse white dwarfs. 

Two EOS have been used for the collapsed objects. One 
of them is one of the standard cold neutron-star EOS. The 
other is a hot EOS suitable for protoneutron stars, which 
are characterized by their high temperature and high lepton 
fraction. 

The equilibrium structure of the collapsed objects in 



all of our models consist of a high density central core of 
size about 20 km, surrounded by a massive accretion torus 
extending over 1000 km from the rotation axis. More than 
90 per cent of the stellar mass is contained in the core and 
core-torus transition region, which is within about 100 km 
from the rotation axis (see Figure [ll]). The central densities 
of the hot protoneutron stars are in the sub-nuclear density 
regime (4 x 10 11 g cm~ 3 <, p <, 2 x 10 14 g cm -3 ). The 
structures of these protoneutron stars are very different from 
those of the cold neutron stars, which the protoneutron stars 
will evolve to in roughly 20 s. The protoneutron stars have 
lower central densities, rotate less rapidly, and have smaller 
values of f3. On the other hand, the structures of the three 
cold neutron stars are similar. Their central densities are 
around 3.5 x 10 14 g cm -3 and their central cores are nearly 
rigidly rotating with periods of about 1.4 ms, slightly less 
than the fastest observed milli second pulsar (1.56 ms). 

Zwerger and Miiller (1997) performed 2D simulations of 



the core collapse of massive stars. The major difference be- 
tween their models and ours is that they used rather simpli- 
fied EOS for both the pre-collapse and the collapsed models. 
When compared with their fastest rigidly rotating model, 
A1B3, we found their pre-collapse star has less total angu- 
lar momentum and smaller f3 than the pre-collapse white 
dwarf of our Model I, although both have the same central 
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density. The differences between their final collapsed models 
(A1B3G1-A1B3G5) and ours are even more significant. The 
values of (3 of our collapsed objects are much larger than 
theirs, suggesting that the EOS plays an important role in 
the equilibrium configurations of both the pre-collapse white 
dwarfs and the resulting collapsed stars. 

The values of f3 of the cold neutron stars are only 
slightly less than the traditional critical value of dynami- 
cal instability, 0.27, frequently quoted in the literature. The 
cold neutron stars may still be dynamically unstable and 
a detailed study is required to settle the issue. Even if they 
are dynamically stable, they are still subject to various kinds 
of secular instabilities. A rough estimate suggests that the 
gravitational-wave driven bar-mode instability dominates. 
The timescale of this instability is about 0.1 s. 
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